A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

11.01.2; 11.02.2 Mkn 421; 11.10.1; 13.07.3; 13.25.3 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



Variability in the synchrotron self-Compton model of 
blazar emission 

A. Mastichiadis, J.G. Kirk 

Max-Planck-Institut ffir Kernphysik, Postfach 10 39 80, D-69029 Heidelberg, Germany 
Received . . . Accepted . . . 



Abstract. We present a model of the spectra of gamma- 
ray emitting blazars in which a single homogeneous emis- 
sion region both emits synchrotron photons directly and 
scatters them to high (gamma-ray) energy before emis- 
sion (a 'synchrotron self-Compton' or SSC model). In con- 
trast to previous work, we follow the full time dependent 
evolution of the electron and photon spectra, assuming 
a power-law form of the electron injection and examine 
the predictions of the model with regard to variability of 
the source. We apply these computations to the object 
Mkn 421, which displayed rapid variability in its X-ray 
and TeV emission during a multiwavelcngth campaign in 
1994. This observation strongly implies that the same pop- 
ulation of electrons produces the radiation in both energy 
bands. By fitting first the observed quiescent spectrum 
over all 18 orders of magnitude in frequency, we show that 
the time dependence of the keV/TeV flare could have been 
the result of a sudden increase in the maximum energy of 
the injected electrons. We show also that different types 
of flare may occur in this object and others, and that the 
energy band most sensitive to the properties of the accel- 
eration mechanism is the X-ray band. 
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1. Introduction 

Mkn 421 is a nearby (z=0.03) BL Lac object, with emis- 
sion characterised by variability at all wavelengths. It ap- 
pears to emit purely non-thermal continuum radiation; 
emission lines are absent. In addition, Mkn 421 is the first 
extragalactic source which has confirmed detections in 
TeV 7-rays (Punch et al. |1992| , Schubncll et al. |1996| ). It 
is also one of the ~50 active galactic nuclei (AGN) to have 
been detected by the EGRET experiment on the Compton 
Gamma- Ray Observatory (von Montigny et al. 1995 ). 



Concerted efforts have been undertaken to obtain si- 
multaneous observations of blazars over a wide spectral 
range. In April and May 1994, Mkn 421 was the target 



of such an observation (Macomb et al. 1995, 1996). This 
object had been observed in previous multi-wavelength 



campaigns (Brodie, Bowyer & Tennant 1987, Makino et 
al. 1987) but the 1994 campaign was the first to cover 



both the GeV and TeV gamma-rays as well as the radio to 
X-ray bands. Furthermore these coordinated observations 
were fortunate to coincide with an increase in the flux 
above 250 GeV by a factor of several in a few days. Con- 
temporaneous observations performed by ASCA showed 
the X-ray flux was also in a high state compared to earlier 
observations. Significantly, EGRET did not detect a flare 
in the MeV-GeV regime, nor was any change observed in 
the radio to optical bands. These results imply that a sin- 
gle population of electrons is responsible for the radiation 
in both the keV and TeV bands and provide a potential 
test for the theoretical models that have been put forward 
to explain the emissi on from AGN (e.g., Maraschi, Ghis- 
ellini & Celotti |l992j, D ermer & Schlickeiser |1993[ S ikora, 
Begelman & Rees 1994 , B landf ord & Levinson 1995 , Mar- 
cowith, Henri & Pelletier $99$) . 



One of the most promising models for the emission 
mechanism of BL Lac objects (including Mkn 421) is the 
"synchrotron self-Compton model" (SSC). Originally ap- 
plied to nonthermal emission of frequency up to X-rays 
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(Jones, O'Dell & Stein 1974), this mechanism has been 
revised to accommodate the recent discoveries by EGRET 
and WHIPPLE. The generic picture is that the radio to 
X-ray photons arise as the synchrotron radiation of rela- 
tivistic electrons accelerated by some unspecified mecha- 
nism in a blob of plasma which itself moves relativistically 
outwards from the core of the AGN; the gamma-rays are 
the result of the inverse Compton scattering of the syn- 
chrotron photons by the relativistic electrons in the blob. 
Several expositions of the steady state spectrum predicted 
by this model can be found in the literature. Variability 
is also predicted in these models, where it arises not from 
a non-steady electron distribution but from changes in 
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the parameters of the source as it moves through the jet 
(Celott i, Ma raschi & Treves |1992], Maraschi, Ghisellini & 
Celotti 



1992, Bloom & Marscher 



1996) 



However, during a flare, the electrons may not have 
time to reach a steady state since their cooling time can 
be much longer than the timescale of the flare. This in turn 
means that in order to follow the evolution of a multifre- 
quency flare such as detected during the 1994 campaign, 
time dependent calculations, taking into account electron 
injection and cooling and following the evolution of the 
spectrum, are indispensable. In this paper we present such 
a calculation. Section 2 contains a brief description of the 
computational method; in Sect. 3 we compare the predic- 
tions of our model to the multiwavelength spectrum of 
Mkn 421 and show the time dependent evolution of flares 
obtained by varying certain parameters in the injected 
electron spectrum. We conclude in Sect. 4 with a discus- 
sion of our results. 

2. Computational model 

As part of a general approach to the modelling of AGN 
spectra, including the self-consistent acceleration of a 
hadronic component, we have developed a numerical code 
capable of following the time evolution not only of the pro- 
ton, but also of the photon and elect ron distributions in 
the source (Mastichiadis & Kirk |l995i henceforth MK95). 
To accomplish this, we assumed each component can be 
represented by a spatially averaged distribution, treat- 
ing escape from the source by a simple 'catastrophic loss' 
term. The system of three time-dependent kinetic equa- 
tions for the distributions as functions of energy is then 
amenable to numerical solution. The relevant physical pro- 
cesses taken into account for the electron and photon dis- 
tributions include synchrotron radiation, inverse Comp- 
ton scattering both in the Thomson and Klein Nishina 
regimes, as well as photon-photon pair production and 
synchrotron self-absorption. Other processes included in 
the original code, such as pair annihilation and photon 
downscattering on cold electrons, turn out to be irrele- 
vant for the case we are considering here. 

Two minor modifications must be made before this 
code can describe the SSC model for blazars. Firstly, since 
the acceleration mechanism is not directly addressed, it is 
not necessary to follow the proton distribution. Instead, 
an arbitrary dependence of the electron injection function 
on energy and time can be implemented. Secondly, the 
extremely high flux of gamma-rays observed from blazars 
indicates that the source must be Doppler-boosted. Thus, 
quantities such as the photon spectrum, which are com- 
puted in the rest frame of the source, require a Lorentz 
transformation into the observer's frame. 

With these changes, the parameters which specify the 
source are as follows: 

— the Doppler factor 5 = l/[r(l — /3cos0)], where T and 
c(3 are the Lorentz factor and speed of the blob, and 



9 is the angle between its direction of motion and the 
line of sight to the observer, 

— the radius R of the source (in its rest frame, in which 
it is assumed spherical) from which the crossing time 
t CT — R/c can be denned; the variation timescale in 
the galaxy frame is then given by t val — Rj (Sc), 

— the magnetic field strength B, specified in units of 
the critical magnetic field: b = 5/(4.414 x 10 13 G). 
When combined with R, this determines the mag- 
netic "compactness" which can be defined in anal- 
ogy with the photon compactnesss (MK95) as 1-q = 
<jtRB 2 / (8irm e c 2 ) , where ctt is the Thomson cross sec- 
tion and m the electron rest mass, 

— the electron injection spectrum, for which we take 
Qc = <7o7~ s e -7 / 7max where 7 is the electron Lorentz 
factor. The three scalar parameters used to specify this 
distribution in the following are s, 7 max and the elec- 
tron injection compactness £ e = ^m e c<7T-R 2 d7(7 — 
1)Q C (see MK95) 

— the effective escape time i e sc of relativistic electrons, 
which can be identified as the timescale over which 
adiabatic expansion losses limit the accumulation of 
relativistic electrons within the source. 

In terms of these parameters, the equation governing 
the evolution of the electron distribution n e is: 



dn^, t) , n c{l,t) 



01 



to 



Q e (n e , n^, 7, t) + L c (n c , n~, 7, tjjl) 



Here L c denotes the electron loss terms (i.e., synchrotron 
losses and inverse Compton scattering), while Q e is the 
injection term. In all cases discussed in this paper the 
contribution of photon-photon pair production to both 
photon absorption and electron injection is negligible, so 
that Q e is essentially a function of only 7 and t. The nor- 
malisation used is one in which time is measured in units 
of the light crossing time of the source (of size R) and 
the particle density refers to the number contained in a 
volume element of size a^R. All quantities are measured 
in the rest frame of the blob. The corresponding equation 
for the photon distribution n 7 reads: 



dn 7 (x, t) 
di 



n 7 (x, t) 



; 7 ^n 7 , n c , X, t) + L 7 (n 7 ,n e ,7,$) 



where here Q-y(n-y, n c , x, t) represents the source terms for 
photons of frequency v = xm c c 2 jh (i.e., synchrotron radi- 
ation and inverse Compton scattering) and L 7 (n 7 , n c , 7, t) 
the loss term, which arises from photon-photon pair pro- 
duction and is negligibly small in the current application. 
The source and loss terms are discussed in detail in MK95. 
Using the above parameters in the code it is possible to 
fit to the multiwavelength spectrum as reported in Ma- 
comb et al. ( [1995 , 



1996) during both the quiescent and 



the flaring stages. 
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Fig. 1. The predicted and observed emission from Mkn 421 
during its quiescent state. Squares indicate observations dur- 
ing the quiescent phase. Tri angle s are observations during the 
flaring state (Macomb et al. 1996). For the adopted parameters 
see text. 



3. Fits to the multiwavelength spectrum of Mkn 
421 

3.1. The quiescent state 

The seven parameters £ e , s, 7 max , $, R, b, and t csc are 
strongly constrained by the observations. In the case of 
Mkn 421, both the X-ray emission and the TeV emission 



are rapidly variable (see Macomb et al. 1995, Takahashi et 
al. 1995, Schubnell et al. 1996). Radio frequency photons, 



however, exhibit a quiescent spectrum. Since the time re- 
quired for electrons to cool by synchrotron emission to en- 
ergies at which they emit substantial numbers of infra-red 
photons is excessively long, we assume that either adi- 
abatic expansion of the blob or escape from the source 
intervenes to limit the time over which a given electron 
radiates in this frequency range. Then the observed infra- 
red to radio spectral index a « 0.35 is directly related to 
s by the well-known formula a — (s — l)/2. 

Although variable, the general form of both the X-ray 
and gamma-ray spectra indicates that in the SSC model 
synchrotron photons are emitted up to a maximum fre- 
quency of approximately 10 18 Hz while inverse Compton 
photons are present up to at least 10 27 Hz. Denoting these 
frequencies by u Sj is x 10 18 Hz and v c ,27 x 10 27 Hz, respec- 
tively, we have 



: « 10-%i 8 

3 X 10V,27 



(3) 
(4) 



(MK95). Eq. (4) assumes that photons of frequency f c ,27 
are produced by inverse Compton scatterings in the Klcin- 



Nishina regime. From these expressions we can immedi- 
ately deduce 7 max and the magnetic field strength in terms 
of the (as yet undetermined) Doppler factor 5: 



7max = 3 X 10 6 ^ C: 27^ 



-1 



B = 5 x 10 3 <5i4,i8f c 27 gauss 



(5) 
(6) 



The next step is to relate the observed bolometric non- 
thermal flux F of the source to the power injected in rel- 
ativistic electrons. Using the normalisation of MK95 we 
require the injected electron compactness £ c in the rest 
frame of the source to be 



3a T FDl 

¥r 



m c c° 



(7) 



where is the luminosity distance to the source given 
by £> L = 2c[z + 1 - (z + l) 1/2 }/H for a q = 1/2, 
A = cosmology; for the Hubble constant we assumed 
H = 75 h km s" 1 Mpc -1 and the redshift of Mkn 421 is 
z — 0.0308. Observations of Mkn 421 indicate that in the 
SSC model, the flux of synchrotron photons (i.e., those 
of frequency < j^is) is comparable to the flux of inverse 
Compton photons (those with frequency > v Sl m). In the 
rest frame of the blob, this implies approximate equality 
of the two compactnesses £ c and £b since these quanti- 
ties are proportional to the internal photon and magnetic 
energy density respectively. Writing rj = i e / is we have 



R = 3 x IQ^l^B' 2 cm 



(8) 



with B in gauss. Using Eq. (7) and the observationally 
derived relation FD^ « 6 x 10 44 erg/sec, we obtain from 
Eq. (8) 



R~1.2x 10 



18 



FDl 



6 x 10 44 erg/sec 



1/2 



V 



-l/2 B -l 5 -2 cm(9) 



The observation that the synchrotron spectrum steep- 
ens at a frequency somewhere between millimeter and op- 
tical wavelengths (Brodie, Bowyer & Tcnnant 1987) en- 



ables one to estimate the escape time first by calculating 
an effective Lorentz factor 7b r below which escape is faster 
than cooling 



Tbr 



(10) 



and then relating this to the turnover frequency by 

1.3 x 10 20 <5&7 2 r . (11) 



This approach gives an escape time 



tcsc = 8.9 X 10 14 5 1 /2 B -3/2 i? -l 1/ -l/2 3 



(12) 



expressed in units of t CI . Here ^,15.3 — fb/10 15 ' 3 Hz is the 
frequency of the spectral break. 
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Fig. 2. Evolution of the multiwavelength spectrum of Mkn 
421 in the case where the amplitude of the electron injection 
changes impulsively by a factor of 3. The solid line corresponds 
to the quiescent spectrum of Fig. 1. The short and long dashed 
lines show the spectrum at 1 and 2t var after the change in the 
electron injection (corresponding to 1.2 and 2.4 days respec- 
tively). The dotted line shows the new steady state which, 
however, is achieved many i var later. 



Finally, the rapid variability of the X-ray and TeV 
emission on timescales i var of 10 5 sees constrains the over- 
all size of the source and the Doppler boosting factor. Thus 
using i var = R/(cS) and eliminating B between Eqs. (6) 
and (9) to find the value of R, we deduce 



15* 



-1/4 



(13) 



where t& — i V ar/(10 5 s) and we have inserted canonical 
values for the observed quantities. 

Following the above guidelines one can allow the code 
to reach a stationary state and find detailed fits to the 'qui- 
escent' multiwavelength spectrum of Mkn 421 (Macomb et 
al. 1995, 1996). With the exception of the slope of injected 
electrons s, all other parameters can be expressed in terms 
of the Doppler factor S. Thus, we first use the approximate 
observed quantities 14, is, ^c,27, FD^, a, fb,i5.3, and the 
observed approximate equality of £ e and £b, to fit the ob- 
served quiescent spectrum and subsequently adjust 6 to 
fit the observed variation time t$. We find a satisfactory 
fit using the parameters R = 4.7 lO 16 ^ 3 cm, B = .07<5i5 
gauss, 7 max = 2 lO 5 ^ 1 , s = 1.7, £ c = 1.93 lO -5 ^ 1 and 
t csc = 3.3t cr (5i 5 . Here S15 = 5/15. Fig. |l| shows the calcu- 
lated quiescent spectrum for <5i5 = 1. 




Fig. 3. Plot of the flux at various frequencies (normalised to 
its quiescent value) for the flare that corresponds to a change 
in q e by a factor of 3 over the quiescent state. The dotted line 
corresponds to a wavelength of 0.4 cm, the small dash line to 
optical wavelengths, the large dash line to 2-10 keV X-rays, 
the small dot-dash line to .1-30 GeV 7-rays while the large 
dot-dash line to 1 500 GeV 7-rays. 



3.2. The flaring state 

In the present model, the simplest way to explain th e flar- 
ing activity of Mkn 421 reported by Macomb et al. ( 1995 ) 
is to change the electron parameters q c and/or 7 max or 
to change the magnetic field. The first two can be under- 
stood as sudden variations in the acceleration mechanism 
while the third corresponds to the blob entering a region 
of higher field strength. 



3.2.1. Flares due to changes in q c 

Fig. U shows how a flare evolves in the case in which the 
amplitude of the electron injection spectrum q e is suddenly 
increased by a factor of 3 above the quiescent value found 
in the previous section and then left constant. The short 
and long dashed lines correspond to snapshots of the spec- 
trum after 1 and 2t var , respectively, measured from the 
time at which q c changed (t val ~ 1.2 days). The dotted 
line corresponds to the new steady state, which is achieved 
many days later. It is evident that a change in q c produces 
prompt flaring activity in the UV/X-ray and TeV 7— ray 
bands. This can be attributed to the fact that high energy 
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Fig. 4. Evolution of the multiwavelength spectrum of Mkn 421 
in the case where the maximum energy of the electrons 7 max 
changes impulsively by a factor of 5. The solid line corresponds 
to the quiescent spectrum of Fig. 1. The short and long dashed 
lines show the spectrum at 1 and 2 i var after the change in the 
electron injection (corresponding to 1.2 and 2.4 days respec- 
tively). The dotted line shows the new steady state. 



electrons have shorter cooling times than low energy ones. 
Therefore the newly injected high amplitude electron dis- 
tribution cools first at the high frequency end, producing 
synchrotron X-rays and inverse Compton TeV gamma- 
rays. The lower frequency photons (optical synchrotron 
and GeV gamma-rays) lag behind somewhat. This can be 
better seen in Fig. ||, which shows the lightcurves of vari- 
ous frequency bands as a function of time (the travel cross- 
ing time across the emitting region has not been taken 
into account). The X-rays have the fastest response while 
the optical shows a slower response. Similarly TeV radia- 
tion has a faster response than GeV radiation. Note also 
that the amount by which the synchro-Compton compo- 
nent increases is larger than that by which the synchrotron 
component rises, since the latter varies approximately as 
proportionaLto q% whereas the former is roughly propor- 
tional to q c . u However, this holds only if the change in q e 
lasts long enough for a new quasi steady state to be es- 
tablished. If, for example, the electron injection is turned 
off after one crossing time, this effect will not be observed. 
Finally, it can be seen from Fig. g, that although the pre- 
dicted flare matches the observed increase in amplitude as 
observed by Whipple, it does not reproduce the hard X- 
ray spectral shape. In addition, the model underestimates 
the high X-ray flux by at least one order of magnitude. 



1 In the present model self-consistent electron cooling is in- 
cluded by both synchrotron and inverse Compton losses. This 
breaks the simple linear dependence of synchrotron flux on q e . 
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Fig. 5. Plot of the flux at various frequencies (normalised to 
its quiescent value) for the flare that corresponds to a change 
in 7 m ax by a factor of 5. The dotted line corresponds to a wave- 
length of 0.4 cm, the small dash line to optical wavelengths, 
the large dash line to 2-10 keV X-rays, the small dot-dash line 
to .1-30 GeV 7-rays while the large dot-dash line to I 500 GeV 
7-rays. 



3.2.2. Flares due to changes in 7„ lax 

The evolution of a flare which corresponds to an impul- 
sive increase of 7 max by a factor of 5 is shown in Fig. || 
Once again the short and long dashed lines correspond to 
snapshots of the spectrum after 1 and 2 t cr respectively as 
measured from the time of the sudden change in 7 max - The 
final steady state is shown as a dotted line and it is reached 
after roughly 3i cr . In contrast to the previous case, the pre- 
dicted flaring activity due to an increase in 7 max gives a 
good fit both in X-rays and TeV 7— rays. Figure || displays 
the light curves predicted in frequencies from infra-red to 
hard gamma-rays. The evolution of the flare is very differ- 
ent from the previous case as here the source exhibits an 
outburst only in X-rays and TeV 7— rays while the other 
frequencies remain practically unaffecte d. Th is behaviour 
was also found by Marscher fc Tr avis ( 1996 ) and it was 
suggested by Macomb et al. ( 1995 ) as a possible interpre- 
tation of their data. Note that the variability in the GeV 
range is of the same order as the size of the error bars. 
Furthermore we find that, contrary to the previous case, 
the X-ray flare is of higher amplitude than the TeV one. 
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Fig. 6. Plot of the 2-10 keV photon index as a function of the 
flux at the same energy band that corresponds to a change in 
7max by a factor of 5. The dotted line corresponds to an impul- 
sive change in 7 ma x, while the full and dashed lines correspond 
to a more gradual change with the change completed in the 
first case in one t wal and in 2t va r in the second. The triangles 
and squares indicate the values at 1 and 2t va _ r respectively. 
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Fig. 7. Plot of the ^,500 GeV photon index as a function of the 
flux at the same energy band that corresponds to a change in 
7max by a factor of 5. The dotted line corresponds to an impul- 
sive change in 7 max , while the full and dashed lines correspond 
to a more gradual change with the change completed in the 
first case in one i var and in 2t var in the second. The triangles 
and squares indicate the values at 1 and 2t var respectively. 



From Fig. [| it is also evident that in this case the out- 
burst is strongest in the X-ray regime. In fact, there are 
marked changes in the spectral index of the X-rays dur- 
ing a flare. This is displayed in Fig. |(| which plots the 2- 
10 keV spectral index as a function of flux in this energy 
range. Three different flares are plotted, which differ in 
the time scale over which 7 max is increased: first of all im- 
pulsively (dotted line), then on a time scale of t vax (solid 
line) and then 2i var (dashed line). In each case a spec- 
tral hardening in X-rays with increasing flux is predicted, 
in qualitative agreement with recent ASCA observations 
(Takahashi et al. 1995| ). A similar effect can be observed 
in the hard gamma-ray (> 500 GeV) range. This is shown 
in Fig. 0. 

3.2.3. Flares due to changes in B 

As a final example, we present a flare caused by a sudden 
increase in the strength of the magnetic field. Fig. |^ de- 
picts the spectral evolution that corresponds to an impul- 
sive increase of the magnetic field strength by a factor of 
3 and then left there. In this case, the whole synchrotron 
spectrum is shifted to the right by the same factor and 
is also boosted in intensity, as might be expected. The 
synchro-Compton component, on the other hand, is re- 



duced compared to its original steady state value. The 
reason is that the ratio £ e /@B is reduced when the field 
increases, so that the total luminosity (which is held con- 
stant) is redistributed towards the synchrotron compo- 
nent. This behaviour can also be seen in Fig. ^| where 
one can observe a fast increase in the flux of low fre- 
quency bands along with a decrease of the flux in the 7-ray 
bands. This result contrasts strongly with similar investi- 
gations by Bloom & Marscher ( 1996| ). There is, however, 
no discrepancy, since these authors consider a change in 
the magnetic field whilst holding the electron distribution 
constant. In our case, we keep the electron injection (and 
hence the total luminosity) constant. 

4. Discussion 

In the present paper we obtained fits to the X-ray/TeV 
flares observed during the 1994 multiwavelength campaign 
of Mkn 421 in the context of the homogeneous synchro- 
Compton model. This model assumes that the most im- 
portant photon targets for inverse Compton scattering 
by relativistic electrons are the synchrotron photons they 
themselves produce. In the case of BL Lac objects such 
as Mkn 421, this assumption may be justified, since there 
is no evidence of a more important photon component. 



A. Mastichiadis, J.G. Kirk: Variability in the synchrotron self-Compton model of blazar emission 



7 




15 20 25 

log V (Hz) 



Fig. 8. Evolution of the multiwavelength spectrum of Mkn 421 
in the case where the magnetic field strength changes impul- 
sively by a factor of 3 and then left constant. The solid line 
corresponds to the quiescent spectrum of Fig. 1. The short and 
long dashed lines show the spectrum at 1 and 2 t var after the 
change in the electron injection (corresponding to 1.2 and 2.4 
days respectively). The dotted line shows the new steady state. 



However, this may not be the case in other sources, where 
photons from an accretion disk (Dermer, Schlickeiser & 



Mastichiadis 1992) or photons sca ttered from broad line 
clouds (Sikora, Begelman & Rees 1994 ) may dominate. 
Nevertheless, a combination of synchrotron radiation and 
inverse Compton scatte ring o n external photons (Dermer, 
Sturner & Schlickeiser 1996 ) may be useful in modelling 
results such as those of the multi-wavelength campaign 
on 3C273 (Lichti et al. 1995). Whether a homogeneous 



synchro-Compton model such as presented here can ex- 
plain these or simila r observations (e.g., those of 3C279: 
Maraschi et al. 1994 Hartman et al. 1996 ) is currently un- 
der investigation. 

We obtain the full time dependent behaviour of flares 
by fitting the 'quiescent' spectrum of the source and then 
varying one of the free parameters. Three simple ways of 
producing flares were investigated: (i) by changing the am- 
plitude of the electron injection spectrum, (ii) by chang- 
ing the maximum energy of the electron injection spec- 
trum and (iii) by changing the magnetic field strength. 
We found a good fit to the observations using a flare of 
type (ii). This produces changes only in the X-ray and 
TeV bands, leaving all the other bands essentially unaf- 
fected. It also reproduces qualitatively the observed hard- 
ening of the X-rays with increasing intensity (Takahashi 
et al. 



1995|) 



We also found that X-rays are first to re- 
act to any change in the electron injection. This is par- 
ticularly pronounced for flares of type (i) and (ii) (see 
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Fig. 9. Plot of the flux at various frequencies (normalised to 
its quiescent value) for the flare that corresponds to a change 
in magnetic field strength by a factor of 3. The dotted line 
corresponds to a wavelength of 0.4 cm, the small dash line to 
optical wavelengths, the large dash line to 2-10 keV X-rays, 
the small dot-dash line to .1-30 GeV 7-rays while the large 
dot-dash line to 1 500 GeV 7-rays. 

Figs. H and |J). A change in the acceleration parameters is 
tracked more closely by X-ray photons than by photons 
in other wave-bands, since the X-ray producing electrons 
have the highest energy and, therefore, the the fastest 
cooling timescale. 

The parameters we find for the fits are similar to 
those found by other authors (see, for example, Marscher 
& Travis fll996j ), Sambruna, Maraschi & Urry ( |1996| )). 
The fast time variability of Mkn 421 (~1 day) implies a 
Doppler factor of ~15. However, a faster variation can eas- 
ily be accommodated in our model, since new fits with a 
larger Doppler factor can be obtained simply scaling the 
parameters as indicated in Sect. 3. One potentially inde- 
pendent constraint on the model is provided by the effects 
of synchrotron self-absorption. At the lower radio frequen- 
cies, the absorption turnover can be seen in Figs. [I], ||, ^ 
and [| For higher Doppler factors, this effect should dis- 
appear, since the required luminosity is then provided by 
a lower value of £ c . As a result, the electron column den- 
sity of the source is reduced. However, the homogeneous 
SSC model presented here is so compact that the effects 
of induced Compton scattering can also be expected to 
manifest themselves in the radio range (Coppi, Bland- 
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ford & Rees 1993, Sincell & Krolik 1994). A simple es 



timate of the importance of this effect can be obtained 
by evaluating the parameter 7i n d = N c aTRTb T /(m e c 2 ), 
where Xb r is the brightness temperature of the radiation 
in energy units and N e is the total electron density in 
the source. Even if we assume that no thermal electrons 
are present, this parameter exceeds unity for frequencies 
less than roughly 700 MHz, given the parameters of our 
quiescent model. Thermal electrons, however, accumulate 
within the source, although we have not calculated their 
density explicitly, consequently, a modification of the sim- 
ple synchrotron spectrum at gigaherz frequencies may be 
possible. 

Our calculations predict the time-dependent flaring ac- 
tivity to be expected given certain idealised fluctuations 
in the injection spectrum of high energy electrons into a 
relativistically moving blob. Although we do not address a 
specific acceleration mechanism, it is possible to interpret 
the overall picture as one in which a shock front rapidly 
accelerates electrons and leaves a trail of them (i.e., injects 
them) in the plasma streaming away from the shock on its 
downstream side. In our model, we assume the typical di- 
mension of the radiating part of the downstream plasma in 
its rest frame is R. The photon escape time is i cr , which we 
find to be roughly one third of the time taken for electrons 
to cross the emitting region. Thus, assuming electrons are 
carried along by the plasma, our picture would indicate 
rapid (~ c/3) movement of the downstream plasma away 
from the shock front. In a more realistic model, the spec- 
trum of injected electrons should also be calculated in a 
time-dependent manner. In fact, the value of 7 max itself 
should probably be determined by a balance between ac- 
celeration and losses in a self-consistent model. This is 
possible in a hadronic model (e.g., MK95), however, self- 
saturation of accelerated protons requires a very high pho- 
ton density, which would render the source opaque to TeV 
gamma-rays. On the other hand, it is possible to imagine a 
model in which protons saturate at extremely high energy 
(e.g., Mannheim 1993), and inject electrons which produce 
the entire emitted spectrum by the synchrotron mecha- 
nism. The time-dependent properties of flares from such 
a model would, however, differ strongly from those found 
above. A detailed investigation of electron acceleration in 
the presence of losses (Heavens & Meisenheimer 1987) has 



so far been performed only in the steady state, and with- 
out considering either bulk relativistic motion or inverse 
Compton scattering. 

Another simplification we have introduced comes from 
the fact that we have used a spherical homogeneous model 
(for a discussion of the problems inhe rent w ith this model 
see, for example, Bloom & Marscher 1996 ). However this 
allows us to understand better the significance of the var- 
ious physical quantities we are using and, at least for the 
case of Mkn 421, it proved adequate to fit the multiwave- 
length spectrum. On the other hand, the inhomogeneous 
models might be superior to the homogeneous ones in the 



sense that they give better overall fits to the AGN spec- 
tra, however they introduce a number of extra parameters 
making a simple understanding of the results difficult. 
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